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ABSTRACT 

We have used the Spitzer sateUite to monitor the mid-IR evolution of SN 1987A over a 5 year 
period spanning the epochs between days ^ 6000 and 8000 since the explosion. The supernova (SN) 
has evolved into a supernova remnant (SNR) and its radiative output is dominated by the interaction 
of the SN blast wave with the pre-existing equatorial ring (ER) . The mid-IR spectrum is dominated 
by emission from ~ 180 K silicate dust, coUisionally-heated by the hot X-ray emitting gas with a 
temperature and density of ~ 5 x 10^ K and ^ 3 x 10'* cm"'^, respectively. The mass of the radiating 
dust is ~ 1.2 X 10~^ Mq on day 7554, and scales linearly with IR flux. Comparison of the IR data 
with the soft X-ray flux derived from Chandra observations shows that the IR-to-Xray flux ratio, 
IRX, is roughly constant with a value of 2.5. Gas-grain collisions therefore dominate the cooling of 
the shocked gas. The constancy of IRX is most consistent with the scenario that very little grain 
processing or gas cooling have occurred throughout this epoch. The shape of the dust spectrum 
remained unchanged during the observations while the total flux increased by a factor of 5 with 
a time dependence of i'0-87±o.20^ ^/ being the time since the first encounter between the blast wave 
and the ER. These observations are consistent with the transitioning of the blast wave from free 
expansion to a Sedov phase as it propagates into the main body of the ER, as also suggested by 
X-ray observations. The constant spectral shape of the IR emission provides strong constraints on 
the density and temperature of the shocked gas in which the interaction takes place. Silicate grains, 
with radii of ~ 0.2 fim and temperature of T ^ 180 K, best fit the spectral and temporal evolution 
of the ~ 8 — 30 fim data. The IR spectra also shows the presence of a secondary population of very 
small, hot (T > 350 K), featureless dust. If these grains spatially coexist with the silicates, then 
they must have shorter lifetimes. The data show slightly different rates of increase of their respective 
fluxes, lending some support to this hypothesis. However, the origin of this emission component and 
the exact nature of its relation to the silicate emission is still a major unsolved puzzle. 

Subject headings: ISM: supernova remnants - ISM: individual (SNR 1987A) ~ 
ISM: interstellar dust - Infrared: general - X-rays: general 



1. INTRODUCTION 

About 10 years after its explosion on February 23, 
1987, supernova (SN) 1987A has evolved from a super- 
nova, when its radiative output was dominated by the 
release of radioactive decay energy in the ejecta, into a 
supernova remnant (SNR), when its radiative output be- 
came dominated by the interaction of its blast wave with 
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the inner equatorial ring (ER). The ER is located at a 
distance of about 0.7 lyr from the center of the explo- 
sion, and could have been pr oduced by mass loss fr om 
a single rotating supergiant (jHeger fc Langeil I1998D or 
by a merger event i n a binary system that a lso formed 
the two outer rings orris &: Podsiadlows ki 20139). The 
transition from SN to SNR was observed at wavelengths 
ranging f rom radio to X - rays, and is summarized in Fig- 
ure 18 in IBouchet eTall mm . 

In this paper we report on the continuing evolution of 
the ~ 5 — 30 /im spectrum and the 3.6, 4.5, 5.8, 8.0, 
and 24 /im photometric fluxes from SN 1987A, spanning 
the ~ 5 year period from day ^ 6000 until day ^ 8000 
after the explosion. Initial re ports and analysis of the IR 
evolution were pr esente d by IBouchet et alT (|2004 |2006[) 
and iDwek et al.l (|2008l ) . In §2 we present the IR data 
obtained by the Spitzer satellite. The evolution of the 
IR emission and the dust composition are presented in 
§3. In §4 we derive the plasma conditions from the IR 
observations. The comparison of the IR emission with 
the X-ray emission, and the evolution of their flux ratio 
are discussed in §5. A brief summary of the paper is 
presented in §6. 
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Fig. 1. — Spitzer images of SN 1987A at Day 7975 (3.6 - 8 fim) and Day 7983 (24 - 70 ^m). These reverse grayscale images use logarithmic 
scaling. 
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Fig. 2.— Spitzer images of SN 1987A at Day 7975 (3.6 - 8 fim) and Day 7983 (24 - 70 ^lm) 
2. OBSERVATIONS AND DATA REDUCTION 



Spitzer 's observations using the InfraRed Array Cam- 
era (IRAC), the InfraRed Spectrograph (IRS), and 
Multiband Imaging Photometer for Spitzer (MIPS) were 
conducted annually during the first two years (2004- 
2005), and then roughly every 6 months for the remain- 
der of its cryogenic mission. SN 1987A was also observed 
incidentally with IRAC and MIPS by the SAGE survey 



eixner et al.ll2006l ) and other projects during 2005. 
Figures [T] and [2] illustrate the IRAC data at the final 
observations (~ day 8000). The IRAC data are all of 
good quality, i.e., the data arc not adversely affected by 
any artifacts induced by bright sources in the field or 
in prior observations. Initial tests of generating mosaics 
from the basic calibrated data (BCD) showed no signif- 
icant difference from the post-BCD mosaics. Therefore, 
aperture photometry was performed on the post-BCD 



3 



10" 



IRS flux ^ ♦ 


24- fim 
8 /zm 








— ^ 

♦ 




5-8 /.nn 

4.5 /j.m 

3.6 /.im 


♦ 



0,100 



6000 



6500 



7000 
time (day) 



7500 



8000 




0.010 



0.001 



2.0 




day 7954 
day 779B 
day 7554 
day 7333 
day 7137 
day 8804 



10 



A (^^.m) 



0.100 



^ 1.6 
? 

0.010^ / ^ ^ 1.4 



S 1.2 



0.001 1 J , , . , , , I 1 o i , I , , , I , , , I , , , I , , , I 

6600 6800 7000 7200 7400 
A (fJ-rn) time (day) 

Fig. 3. — Top row: The evolution of the IRAC and MIPS 24 /im luminosities (left panel) and the specific intensity of the low-resolution 
spectra (right panel) of the equatorial ring. Bottom row: Spectra normalized to the one obtained on day 7554 (left panel). Because of 
the constancy in the shape of the spectrum, any of the IRAC or MIPS photometric fluxes can be used as a proxy for the integrated IR 
spectrum from the ER. The right panel depicts the almost constant ratio between the IRS integrated over the 5-30 /im spectrum to the 
24 /im luminosity. 



7600 7800 8000 



mosaics in order to obtain the broad band evolution of 
SN 1987A at 3.6, 4.5, 5.8 and 8 /zm. These measurements 
and their statistical uncertainties are shown in Figure [31 

Figures [Band [2] also illustrate the MIPS 24 and 70 /im 
images of SN 1987A at Day 7983. There is no evident 
point source at 70 fim. As for IRAC, there seemed to 
be no advantage to reprocessing the BCD, thus aperture 
photometry was performed on the post-BCD mosaics. 
The 24 fim flux densities are shown in Figure [3] At Day 
7983, a 70 yum upper limit was determined by adding an 
artificial point source to the MIPS image at the loca- 
tion of the SN. We started with a source fiux that was 
clearly visible above the background confusion, and then 
reduced the flux until the source was no longer detected 
using the FIND routine of the IDLASTRO library. The 
results indicate the S'(70/im) < 0.09 Jy. 

A signiflcant complication with the IRAC and MIPS 
photometry is that the SN is not resolved from the nearby 
Star 2 and Star 3. If the stars are assumed to have mid-IR 
magnitudes that are equal to the K magnitudes (15.06, 
15.80 respectively from Walborn et al 1993), then their 
combined IR flux densities should be 0.41, 0.26, 0.16, 
0.09, and 0.01 mJy at 3.6, 4.5 5.8, 8, and 24 /im, respec- 
tively. These assumptions may be accurate for Star 2 , 
which is identifled B2 III star (iWalborn et al.lfl99l 
iScuderi et an il99d). However IWalborn et all (|1993[ ) re- 
port both, an IR excess and a variability by at least 
Am = 1 mag for Star 3. The fluxes in Figure [3| are 
plotted after subtraction of the emission of Stars 2 and 



3 under these assumptions. 

IRS observations were reduced starting with the BCD, 
and using the basic steps outlined in the IRS Data 
Handbook. For the low resolution (i? ^ 60 — 120) 
data, we constructed a 2-D background frame using 
the observations of a specific off-source target position. 
IRSCLEAN_MASK was used to identify and clean rouge 
pixels in the background-subtracted BCD frames. The 
cleaned BCD data were averaged for each nod (positions 
of the source at 1/3 or 2/3 the length of the slit), and 
spectra were extracted from these super BCDs us ing the 
Spitzcr IRS Custom Extraction (SPICE) softwar(0: Af- 
ter extraction, the spectra orders were merged to produce 
the full spectra shown in Figure [3|(top right panel). The 
high resolution (R ~ 600) data were processed in a simi- 
lar fashion. The major difference was that instead of us- 
ing a background from a dedicated off-source target, the 
background is derived from the average of two pointings 
that flank the SN on opposite sides (15" and 25" away for 
the short high resolution (SH) and long high resolution 
(LH) spectra respectively). The position angle of these 
pointings varied as dictated by the scheduling of the ob- 
servations, but the distance of the pointings from the SN 
remained constant. For both the high and low resolution 
data, we used the Spectroscopic Modeling Analysis and 
Reduction TooO (SMART) to flt and extract line fluxes. 



http: //ssc . spltzer . caltech . edu/dataanalysistools/tools/sp lce/j 
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Fig. 4. — Left panel: Single composition fit to the IR spectrum of the ER. The fit corresponds to the spectrum of astronomical silicate 
grains radiating at 180 K. The excess 5 — 8 ^m emission evolves at the same rate as the silicate emission (see Fig. [Sj strongly suggesting that 
it is generated by an additional dust component in the ER. Right panel: The addition of any dust component with a smooth spectrum 
which will not reduce the goodness-of-fit to the 9.7 and 18 ^m silicate features is a viable candidate for this second dust component. Here 
the spectrum is fit with a modified blackbody with a A^^ cmissivity law and a temperature of 346 K. 



TABLE 1 

Observed X-ray and Infrared Fluxes From SN 1987A-'- 



day2 


X-ray fiux'^ 


IR fiux-i 


IRX^ 


6067 


(1.21 ±0.18) X IQ-^^ 


(3.15 ± 0.46) X 10-12 


2.6 ±0.5 


6184 


(1.40 ±0.19) X 10-12 


(5.38 ± 0.37) X 10-12 


3.8 ±0.6 


6552 


(2.72 ± 1.0) X 10^12 


(7.43 ±0.38) X 10-12 


2.7 ± 1.0 


6734 


(3.47 ±0.9) X 10-12 


(8.52 ± 0.40) X 10-12 


2.4 ± 0.7 


6829 


(3.62 ±0.79) X 10-12 


(9.07 ±0.38) X 10-12 


2.5 ±0.6 


7159 


(4.43 ±0.66) X 10-12 


(1.12 ± 0.04) X 10-11 


2.5 ±0.4 


7310 


(4.84 ±0.90) X 10-12 


(1.22 ± 0.04) X 10-11 


2.5 ±0.5 


7490 


(5.10 ±0.65) X 10-12 


(1.32 ± 0.04) X 10-11 


2.6 ±0.3 


7690 


(5.83 ± 1.18) X 10-12 


(1.43 ±0.04) X 10-11 


2.5 ±0.5 


7983 


(5.79 ± 1.00) X 10-12 


(1.55 ± 0.04) X 10-11 


2.7 ± 0.5 



1 All fluxes are in units of erg cm— 2 s— i. 

2 Measured since the explosion. 

^ Soft X-ray fiux in the 0.5-2.0 keV band, interpolated to the epochs of the Spitzer/MIPS observations and corrected for an extinction 
column density of Nh = 2.35 X 102l cm-2. 

* IR fluxes on days 6184 to 7983 were derived from the MIPS 24 fim data using the conversion factor of 1.64 (see Fig. [3J- The fiux on day 
6067 was derived using the 10.4 ^m T-ReCS fiux as a proxy for the integrated 5-30 fira spectrum. 
^ The ratio of the IR to soft X-ray fiux from the ER. 



3. THE EVOLUTION OF THE INFRARED EMISSION AND 
THE DUST COMPOSITION IN THE ER 

Figure [3] depicts the evolution of the IRAC and MIPS 
broad band photometry (top left panel), and the low- 
resolution IRS spectrum of SNR 1987A (top right panel), 
as a function of time. The relatively strong background 
may affect the continuum at A > 20 /im and the [Ne III] 
15.6 fim and [S III] 18.7, 33.5 fim emission lines. Line 
emission from SN 1987A and its surrounding medium will 
be presented elsewhere (Arendt et al; in preparation). 
The figures show that the IR flux is rising smoothly in 
all bands, increasing by a factor of ^ 3 between days 6190 
and 7980. The bottom left panel of the figure shows the 
spectra normalized by least-square fit to the one obtained 
on day 7554. The spectrum maintained a nearly constant 
shape throughout the 1800 days of observations. This 
constancy enables the use of any IRAC or 24 /im MIPS 
band as a proxy of the total integrated IR flux from the 
ER. The bottom right panel in the figure shows the ratio 
between the integrated IR flux and the 24 fim flux from 
the ER. The ratio is constant with a value of 1.64 to 
within < 10%. The constancy in the spectral shape also 



places strong constraints on the parameters of the hot 
X-ray gas in which the grains are embedded. 

Figure |4] (left panel) shows that the mean spectrum of 
the ER can be well fltted by astronomical silicate grains 
radiating at a single temperature of ~ 180 K. The relative 
ratio between the two peaks of the silicate features pro- 
vide strong constraints on the temperature. The mass of 
the silicate dust is 1.2 x 10^^ Mq on day 7554, and should 
scale linearly with the IR flux for other epochs (see Table 
1). The figure also shows the presence of an IR excess 
at wavelengths between ~ 5 — 8 /xm that cannot be at- 
tributed to emission from the silicate grains. This e xcess 
emission was first detected bv lBouchet et all (|2006l ). but 
was not further pursued because of the shorter integra- 
tion time and larger noise of that early Spitzer spectrum. 
An unidentified K-band (2.2 /im) continuu m at a level of 
^0.4 mJy, detected bv lKjar et al.l ()2007[ ). may be part 
of this emission component. This emission component, 
now clearly present in all epochs, can be approximated 
in the 5 to 8 /im wavelength region by a v"' power law 
with a spectral index oi a ~ —1.8. This rules out ther- 
mal bremsstrahlung, which has a well defined spectral 
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Fig. 5. — Two-composition fits to the observed infrared spectrum of the ER, indicated by the violet curve. An additional hot component 
consisting of either carbon, metallic iron, iron oxide, or iron sulfite grains (red curves) provides a good fit to the 5 — 8 fim segment of the 
spectrum. The blue curve represents the spectrum of the silicate grains, and the green curve represents the sum of the emission from the 
silicate and the secondary dust component. 



index of a « —0.1 in this frequency regime, as a possi- 
ble source. Synchrotron radiation can also be ruled as a 
source of this emissio n component. The ob served radio 
synchrotron emission (jZanardo et all 120101) has a spec- 
tral index of « —0.8, and its extrapolation from a value 
of 300 mJy at 1.4 GHz on day 7000 will make a negligi- 
ble contribution to the 8 /xm emission around that epoch. 
We therefore conclude that this emission arises from an 
additional hot dust component that is tightly correlated 
with the silicate emission. Its spectral shape should be 
featureless at longer wavelengths as to not reduce the 
goodness of the fit to the observed 9.7 and 18 /im silicate 
emission features. The right panel of the figure shows the 
fit of a blackbody spectrum, modified by a emissivity 
law and a temperature of ^ 350 K, to this component. 
A emissivity is consistent with metallic or carbon 
grain co mposition, rather than diel ectric or silicate com- 
position (jBohren fc Huffmanlll983() . 

The modified blackbody is of course an idealized rep- 
resentation of this secondary dust component. Figure [5] 
shows a two-component fit to the IR spectrum when the 
secondary component was taken to consist of either car- 
bon, pure iron, magnetite (FC3O4), or iron sulfide (FcS) 
grains. All components, combined with the main silicate 
dust component, provided acceptable fits to the over- 
all IR spectrum of the ER. The temperature of the sec- 
ondary dust component is significantly higher than that 
of the silicate grains, ranging from ~ 370 K for pure 
iron grains to 460 K for carbon grains. Some other 
possible secondary dust components, in particular iron- 



and aluminum-oxides, were ruled out because these grain 
species produce strong emission features at wavelengths 
longward of ^ 8 /im that are not present in the Spitzer 
IRS spectra. 

The problem of accounting for the widely different tem- 
peratures of the two dust components, even though they 
may be collisionally-heated by the same gas, will be ad- 
dressed in §4 below. 



4. PLASMA CONDITIONS AS DERIVED FROM IR 
OBSERVATIONS 

Figure [S] shows the plasma densities required to colli- 
sionally heat the different grain types to their respective 
temperatures given in Figure [5] as a function of grain 
radius. The grains are assumed to be embedded in an 
ionized plasma with a temperature of 5 x 10^ K, reflecting 
the average temperature of the soft X-ray component (see 
§5 below). If we require both dust components to arise 
from the same gas, then only the carbon grains support a 
range of plasma densities (enclosed between the two hor- 
izontal dashed lines) that overlap with that of the silicate 
grains. So the IR spectrum of the ER is consistent with 
the emission from a population of small (a < 0.05 /xm) 
carbon grains intermixed with a population of large 
(a > 0.20 fj,m) silicate grains embedded in a plasma with 
a density rig « Uh ~ (2 — 4) x 10* cm~'^. Requiring a 
common environment for the two dust components rules 
out the other dust compositions (Fe, Fe304, and FeS) as 
viable dust candidates. Alternatively, the emission could 
arise from these non-carbonaceous grains, provided they 
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reside in a significantly denser phase (rie ^ 2 x 10^ cm~'^) 
of tlie X-ray emitting gas. Tliis pfiase tlien lias to be 
tightly correlated with the lower-density one in order to 
explain the tight correlation between the IR light curves 
of the two dust emission components. 

In principle, the second dust component could consist 
of grains that have formed in the SN ejecta. There is 
ample evidence for t h e formation of du s t in SN1987A 
(iMoselev et al.l 119891: iLucv et al.l 119911: IWooden et all 
I1993D . In particular, the disappearance of IR fine struc- 
ture line emission from Fe has been interpreted as ev- 
idence for the formation of ir on dust in the ejecta 
(jPwek et al.lll992l:IWoodenlll997[) . However, it is unlikely 
that this second dust component is ejecta dust. The to- 
tal IR luminosity of this second dust component is about 
10^^ erg s~^, calculated for a distance of 50 kpc to the SN. 
This luminosity is higher by a factor of 10^ from that 
expected from the radioactive decay of ^^Ti, assuming 
that ^ 10~^ Mp) of " ^^Ca was produced in the explosion 
(jWooslev et al.lll989D . So the second dust component can 
only be heated by a reverse shock expanding through the 
ejecta. As shown below, the IR light curve depends on 
the rate at which the mass of the dust is swept up and its 
rate of destruction. These rates depend on the medium 
morphology and density, and the shock velocity. The 
tight correlation between the light curves is unlikely to 
arise from separate shocks propagating into vastly dif- 
ferent environments such as the ejecta and the ER. It is 
therefore very unlikely that the emission from the sec- 
ond dust component arises from SN-condensed dust in 
the ejecta. 

The potential presence of carbon with the silicate dust 
in the ER is surprising, defying the common paradigm 
that either carbonaceous or silicate dust can form in the 
wind of the progenitor, depending on the C/0 abun- 
dance ratio in the outflow. The observations therefore 
suggest either that CO formation did not exhaust all the 
carbon in the outflow, or that the co-existence of the 
two dust components is a manifestation of a binary ori- 
gin of the ER. The progenitor of SN 1987A could have 
shared a common envelope with a less massive carbon- 



rich star. The interaction between them could have led 
to the formation of the ER as well as the two outer rings 
observed by the Hubble to be in the vicinity of the SN 
(jMorris fc Podsiadlo^ski|[2009l ). 

5. THE EVOLUTION OF THE INFRARED AND X-RAY 
EMISSION AND THEIR RATIO 

A rise in the IR emission is expected as the shock 
sweeps up more dust during its propagation into the 
dusty circumstellar medium. The rate of increase of the 
IR emission depends on the morphology of the medium 
into which the shock is expanding, the expansion rate of 
the shock, and on the possibility that dust may be de- 
stroyed in the postshock gas. The expected increase of 
the IR emission with time is schematically illustrated in 
Figure [71 The left column represents the case in which 
the grain lifetime is sufficiently long that it can be ig- 
nored, that is, Tsput > t' , where t' is the residence time 
of the dust in the hot gas. The right column, for which 
Tsput < i', represents the case in which the grain life- 
time is sufficiently short that an equilibrium is estab- 
lished between the rate of grain destruction and the rate 
at which the dust is replenished by the expanding shock. 
For both cases the figure depicts a shock expanding into a 
1-dimensional protrusion (top) , a 2-D disk, or ring (mid- 
dle), and a homogeneous sphere (bottom). The depen- 
dence of the volume of the swept-up gas is presented for 
two different phases in the evolution of a blast wave: (1) 
the free expansion phase, during which the velocity, Vf, 
is constant; and (2) the Sedov- Taylor phase, which com- 
mences after most of the kinetic energy of the explosion 
has thermalized, during w hich the velocity Vg ^ t'^^^^ 
(iZel'Dovich fc Raized[l967h . 

When Tsput > t', following the left column of Figure [3 
the intensity of the IR emission scales as the mass of 
swept up dust, which for a constant dust-to-gas mass ra- 
tio, simply scales as the volume of the shocked gas. For 
a freely-expanding blast wave, the volume of the shocked 
gas, Vf, will increase as t' , i'^, and t'^, for the three 
morphological scenarios, respectively. During the Sedov 
phase the volume of the shocked gas, Vs, will be increas- 
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Fig. 7. — A schematic view of the evolution of the volume of shocked swept-up dust for different geometries of the circumstellar medium: 
a linear protrusion (top row); a hollow disk (middle row); a uniform sphere (bottom row). The amount of radiating hot dust is depicted 
as a dark-shaded region behind the advancing shock. The left column depicts the evolution when grain destruction is insignificant, i.e. 
when Tsput > t, where t is the time when the shock first encounters the circumstellar medium. The mass of shock-heated dust increases in 
proportion to the swept up gas. The right column depicts the volume of the shock-heated dust when Tsput > t. The amount of radiating 
dust is then constant as the increase in the mass of the swept-up dust is balanced by its destruction. Each scenario exhibits a different 
dependence on time, which also depends on the expansion law of the blast wave; Vf is the volume swept up for a freely-expanding blast 
wave {vs ~ t"), and Vs is the volume swept if the blast wave entered the Sedov phase of its evolution [vs ~ i""^/^). 
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Fig. 8. — The evolution of the Prenormalized) IRAC and MIPS 
24 fiva intensities as a function of r = t — io , where to is the time at 
which the SN blast wave first encountered the ER. The value of to 
is equal to 5700 d for the IRAC fluxes, and equal to 5500 d for the 
24 ^m flux. The total IR flux emitted by the silicate dust (repre- 
sented by the evolution of the 24 fim luminosity) exhibits a steeper 
rise than that of the secondary dust component, represented by 
the IRAC luminosities. The rise in the IR spectrum is consistent 
with a Sedov blast wave expanding into a disk-like medium, with 
negligible grain destruction see Figure [7] 

ing of increasing more slowly as t'^/'^, t''*/^, and t'^/^, 
respectively. 

When Tsput < t' , following the right column of Fig- 
ure [21 the dust column density through the shock will 



eventually reach a constant value, roughly equal to 
Vs X Tsput, where Tsput is the grain destruction time in 
the postshock gas. The IR intensity behind a freely- 
expanding shock will therefore reach a constant value if 
the shock is expanding into a one dimensional protru- 
sion, and will increase as t' or when it expands into a 
ring or sphere, respectively. During the Sedov phase the 
volume of the shocked gas, Vs, will actually be decreas- 
ing as t'""^/^, or increasing more slowly as t'^/^, 
for the three morphological scenarios, respectively. 

Figure [8] depicts the evolution of the IRAC and 24 /im 
MIPS fluxes. The fluxes were fitted by a power-law fit of 
the form (t — io)"j where t is the time since the explosion, 
and a and to are free parameters. The value of to can be 
interpreted as the time when the shock first encountered 
the ring. The values of {a, to} were found to be {0.6, 
5700 d} and {0.87, 5500 d} for the IRAC and MIPS 
fluxes, respectively. 

The dust flux scales as the volume, and the value of a 
derived from the light curve is consistent with Vs ^ t''^!^ 
(left column of Figure [T]) , indicating that grain destruc- 
tion is unimportant. More specifically, the 24 /im flux in- 
creases at a rate that is consistent with an Sedov- Taylor 
blast wave expanding into a circumstellar disk/ring. As 
shown below, the inferred dust lifetime is consistent with 
currently used values for the sputtering rates of dust 
grains in hot plasmas. 

The value of to is consistent with the time inferred from 
X-ray observations that showed a deceleration of the ra- 
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Fig. 9. — Left panel: The evolution of the IR flux, represented by the 24 fj,m MIPS (filled blue diamonds) and 10.4 fim T-ReCS (open 
blue diamond) observations, the soft X-ray flux (derived from the Chandra observations), and IRX, the IR-to- X-ray flux ratio. Infrared 
and X-ray fluxes are given in units of 10"^'^ erg cm~'^. IRX is dimensionless. Right panel: Details of the evolution of IRX. Solid blue 
diamonds and the open green diamond are derived using the MIPS 24 fim and the T-ReCS 10.4 iim fluxes as proxies for the integrated IR 
spectrum. The IRX derived using the integrated IR fluxes is shown as fllled red diamonds. The data are consistent writh a constant value 
of IRX over the entire ~ 6000 - 8000 d epoch. 
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dial e xpansion of the SN around day 6000 (jRacusin et alJ 
l2009f l. and the upturn in th e soft X-ray lumin osity that 
occurred around day 6200 (jPark et al.|[200l . The IR 
and X-ray observations therefore offer a consistent sce- 
nario in which both fluxes originate from the dust and 
gas shocked by the blast wave that has started to pen- 
etrate the main body of the ER around day 6000 after 
the explosion. 

5.1. The Evolution of IRX 

Analysis of the Chandra X-ray observations of the ER 
show that the spectrum is generated by two distinct emis- 
sion components: a soft component with temperatures 
~ (0.3 — 0.6) kcV, and a harder component characterized 
by temperatures betw een ~ (2 — 5) keV (jZhekov et alJ 
[200llPark eraI1[2006l) . The emission from the soft com- 
ponent arises from the gas heated by the shock that is 
transmitted into the ER. The hard component arises 
from the lower density gas interior to the ER that is 
shocked twice: once by the advancing SN blast wave, 
and a second time by the shock reflected from the ER. 
We assume here that the IR emission arises from the soft 
X-ray component generated by the shock that is pene- 
trating the ER. 

Figure O (left panel) depicts the evolution of the soft 
X-rays and IR fluxes, and their ratio, IRX, as a func- 
tion of time. The IRS observations only cover the 
~ 6804 - 8000 d epoch, so we used the MIPS 24 fim 
observations that cover the additional ~ 6190 — 6800 d 
epoch as a proxy for the integrated 5-30 /im emission. 
The figure also depicts the day 6067 datum point, for 
which we used the scaled mid-IR p hotometry from Gem- 
ini T-ReCS 10.4 fiui observations (|Bouchet et all 120061 ) 
as a proxy for the integrated IR flux. 

The radial expansion of the soft X-ray component 
slowed down around day ~ 6000, which was interpreted 
as evidence for the transition of the blast wave from the 
free-expansion to the Se dov- Taylor phase as it enters the 
maiii body of the ER (|Park et al.l 120051 : iRacusin et alJ 
l2009f ). The T-ReCS datum suggests that this transition 
has also manifested itself in the IR emission, re-enforcing 
the consistency between the X-ray and IR evolution of 
the remnant. 



The evolution of IRX is shown in more detail in the 
right panel of the figure. It has a roughly constant value 
of ~ 2.5 throughout the entire ~ 6000-8000 d epoch. An 
exception is the 6190 datum point. We believe that the 
anomalously large value of IRX at this point is caused 
by a relative dip in the X-ray flux starting near day 5800 
and ending with the rise in the soft X-rays near day 6200. 
The large value of IRX, also suggests that a significant 
fraction of the refractory elements in the ER are locked 
up in dust. 

A constant value of IRX indicates that the IR and X- 
ray fluxes evolve at the same rate. This requires the dust- 
to-gas mass ratio to be uniform throughout the medium 
into which the shock is expanding, and that neither grain 
destruction nor gas cooling affect the IR and X-ray emis- 
sion from the shocked gas. Alternatively, the constancy 
in the ratio could be the result of both processes, grain 
destruction and gas cooling, being important. In this 
case, the IR and X-ray emission have attained a steady 
state in which the losses to grain destruction and cool- 
ing are balanced by the influx of new material such that 
the volume of radiating gas and dust remains constant 
in time. However, the rate of increase in the IR fluxes 
suggest that the volume of the shock heated dust is in- 
creasing, so that the second scenario is a-priori not likely 
to be the cause for the constancy in IRX. In the follow- 
ing we examine in more detail the physical condition of 
the shocked gas that can lead to a constant IRX. 

5.2. Grain Destruction in the Hot Gas 

The lifetime of a dust grain of radius a moving at a high 
velocity through a hot g as with t emperatures bet ween ^ 
10^- 10'^ K is given by (jNozawa e t al. 20il6i; ,Dwek et alJ 
I2OOI : 

/ \ A a{fim) 

Tsputiyr) = Asput (1) 



ni{{cm ■^) 



where A 



sput 



3.3xl0^9.0xl0^ 5.6xl0^ 2.8xl0^and 
3.9 X 10^, respectively, for silicate, carbon, iron, FeS, and 
Fe304 dust, calculated for grains moving at velocities 
> 500 km s~^ through a 5 x 10^ K gas. The values of 
the velocity and temperature characterize those of the 
shock giving rise to the soft X-ray component. Figure [TOl 
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Fig. 10. — The sputtering lifetime of the different grain species calculated for the same gas density-grain radii combinations required to 
heat the different dust species to their observed temperature. Dust species are color coded as in Figure [6] Within the permitted range of 
silicate grains, the straight lines indicate the loci of constant sputtering lifetimes Teput = 3, 6, 9 (bold), 12 yr. For carbon dust the lines 
indicate Tsput = 0.4, 0.6, 0.8 (bold), 1.0 yr. For iron oxide dust the hnes indicate Tgput = 0.02, 0.03 (bold), 0.04, 0.05 yr. 




10° 10' 10° 10 20 30 40 

Te (K) time (yr) 

Fig. 11. — Left panel: The equilibrium cooling curve for a plasma with ER abundances (thick line). The thin lines represent the gas 
cooling via dust-gas collisions. All curves assume that all the dust is represented by silicate grains with radii indicated in the figure. Right 
panel: The evolution of the temperature behind shocks of different velocities, expanding through a homogeneous dusty medium with ER 
abundances and a preshock hydrogen density of 1 X 10* cm~^. 



depicts the sputtering lifetime of the dust as a function 
of grain radii and gas density. The figure shows that the 
hfetime of the sihcate grains is between ^ 4 and 15 
yr, depending on the combination of grain size and gas 
density. In contrast, for the secondary dust component 
this same range of gas temperatures and densities yields 
significantly shorter grain lifetimes, between ~ 0.4 and 
~ 1 yr, which are too short to explain the observed tem- 
poral evolution of the secondary dust component. 

5.3. Cooling of a Shock-Heated Dusty Plasma 

Figure [TT] depicts the equilibrium atomic cooling rate, 
A(re), of a gas with ER composition (thick line, left 
panel) as a function of electron t emperature, T p . Th e 
ER abundances were taken from iMattila et al] (|2010( ). 
and were derived from AAT and VLT observations of 
the supernova between days ^ 1400 and 5000 after the 
explosion. The thin lines in the figure show the cooling 
of the plasma via gas-grain collisions, assuming a pure 
silicate dust composition. The silicate dust-to-gas mass 



ratio was assumed to be 0.0029, half the solar value. The 
different curves assume that all the dust is represented 
by grains with the indicated radii. Dust cooling domi- 
nates the cooling via atomic transitions at temperatures 
above ~ 10^ K, depending on grain size. The value of 
IRX is the ratio of these two cooling rates. From the 
fact that the initial value of IRX is only ^ 2.5 we can 
conclude that the grain sizes in the ER must be rela- 
tively large, with typical radii of > 0.2 /im, consistent 
with the sizes required to heat the dust to a temperature 
of ~ 180 K (sec Fig. [6|). Using this lower limit on the 
grain radii in eq. ([T]) , gives a postshock density of about 
(2 — 4) X 10^ cm~^, where we assumed an upper limit of 
0.5 /im on the grain radii. 

The equilibrium cooling may not represent the actual 
cooling of a shocked gas, which may be far from ion- 
ization equilibrium. The right panel in Figure [TT] de- 
picts the temperature profile behind shocks of different 
velocities expanding into a homogeneous dusty gas with 



10 



a prcshock hydrogen density of 10"* cm~'^ and ER com- 
position. The shock models include the effects on non- 
equilibrium ionization (NEI) , and the effect of dust cool- 
ing. The ^ 600 — 700 km s""'^ shocks well represent 
the 3 — 7 keV teniperature of the soft X-r ay emitting 
plasma (iPark et al.ll2"oMlZhekov"eFaI|[2009l ). The figure 
shows that the cooling time of this shocked gas is about 
12 — 20 yr for a preshock gas density of uh = 10^ cm~^. 
The IR observations suggest that the postshock density 
is ^ (2 — 4) X 10"* cm^ , giving a gas cooling time of 
~ 20 - 40 yr. 

5.4. Scenarios for the Constancy of IRX 

The results presented above show that the timescales 
for grain destruction and gas cooling are sufficiently long 
that neither processes has significantly affected the IR or 
X-ray emissions, or their flux ratio, IRX. The sputter- 
ing timescale suggests that grain destruction may become 
important only at day - 6000 -f- 3200 w 9200 after the 
explosion, which is about 25 yr. The cooling time sug- 
gests that the X-ray emission may not be affected until 
t ^ 20 — 40 yr. So the scenario in which both the cooling 
time and the sputtering time are longer than the age of 
the shocked gas, is a viable explanation for the constancy 
of IRX. 

The alternative scenario (see §5.1) suggests that the 
approximate constancy in IRX is the result of both pro- 
cesses, grain destruction and gas cooling, having attained 
a steady-state, in which the amount of prcsliocked gas 
and dust entering the shock is balancing the amount of 
dust that is destroyed and the gas that has cooled. Such 
scenario requires the grain destruction timescale and gas 
cooling time to be significantly shorter than ~ 100 d, re- 
quiring postshock densities that are more than an order 
of magnitude larger than those inferred from the temper- 
ature of the coUisionally-hcated silicate dust. 

The original scenario in which the evolution of the IR 
and X-ray emission from the postshock gas has not yet 
been affected by the two processes offers therefore a more 
reasonable explanation for the constancy of IRX. It re- 
lies on the existence of a population of relatively large, 
a > 0.2 ^m, silicate grains, embedded in a gas with 
postshock densities of about (2 — 4) x 10^ cm~'^. These 
numbers should be regarded as suggestive values for the 
model parameters. 

5.5. The Mystery of the Secondary Grain Component 

However, this simple explanation is complicated by the 
presence of the secondary dust component consisting of 
much smaller dust grains. Figure [10] shows that the 
density-radius combinations required to heat the carbon 
dust to ^ 460 K will also destroy the grains in less than 
half a year, and the various iron bearing grains would 
last for less than ~ 10 d, which definitely rules them out 
as viable candidates for the secondary dust component. 

The short sputtering time for the carbon dust suggests 
that their column density through the shock reached a 
steady-state between injection and destruction shortly 
after the shock penetrated the ER. Consequently, the IR 
emission from these dust particles should have increased 
at a different rate than the emission from the silicate 
grains. This is in contrast to the good correlation be- 
tween the two emission components. 



Figure |5] shows that the 5 — 8 /im emission does increase 
slightly less rapidly than the 24 ^m emission. However, 
the difference in slope is only 1/4 as large as expected 
if the Tsput (hot dust) < t' < Tsp„t (sflicate). It is there- 
fore possible the ~ 5 — 8 yum emission could arise from a 
yet unknown dust component consisting of grains with ei- 
ther radii or IR emissivities (or both) sufficiently small to 
reach the observed high temperature. They should also 
have a featureless spectrum in the '-^4 — 40 /im wave- 
length region, in order not to reduce the goodness of the 
fit to the ~ 9 — 20 fini silicate spectrum. Their bind- 
ing energies should be sufficiently large, or the ambient 
density should be sufficiently low, so that their sputter- 
ing lifetime be longer than 5 yr in spite of their small 
radius. The nature of this secondary dust component is 
still an unresolved aspect of the observations. 

6. SUMMARY 

In this paper we presented the mid-IR evolution of SN 
1987A over a 5 year period spanning the epochs between 
days ^ 6000 and 8000 since the explosion. Its radiative 
output during this epoch is dominated by the interaction 
of the SN blast wave with the pre-existing equatorial ring 
(ER) . The main results of this paper can be briefly sum- 
marized as follows: 

1. The ~ 8 — 30 fiia mid-IR spectrum is dominated by 
emission from ^ 180 K silicate dust, collisionally- 
heatcd by the hot X-ray emitting gas with a tem- 
perature and density of ^ 5 x 10^ K and ~ (2 — 4) x 
10^ cm"'^, respectively. The mass of the radiating 
dust is ~ 1.2 X 10~^ Mq on day 7554, and scales 
linearly with IR flux. 

2. A secondary emission component dominates the 
spectrum in the ^5 — 8 /im region. Its inten- 
sity and spectral shape rule out any possible gas 
or synchrotron emission mechanism as the source 
of this emission. It must therefore attributed to 
a secondary dust component radiating at temper- 
atures above ~ 350 K. 

3. The overall shape of the ~ 5 — 40 /im dust spectrum 
has not changed during the observations, suggest- 
ing that the density and temperature of the soft 
X-ray emitting gas have not significantly changed 
during the more than 5 years of IR observations. 
The constancy in the spectral shape of the IR emis- 
sion also suggests that the mass ratio of the silicate 
to the secondary dust component remained roughly 
constant during this period. 

4. The evolution of the IRAC and MIPS fluxes can 
be described by a power law in time since the 
first shock-ER encounter. The silicate emission 
increases as consistent with X-ray observa- 
tions, suggesting that the blast wave has transi- 
tioned from a free expansion to the Sedov phase, 
and is now expanding into the main body of the 
ER. 

5. The infrared-to- X-ray fiux ratio, IRX, is constant 
with a value of ~ 2.5 throughout this epoch. The 
magnitude of IRX shows that the cooling of the 
shocked gas is dominated by IR emission from the 
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coUisionally- heated dust with radii ^0.2 fim^ and 
that a significant fraction of the refractory elements 
in the ER should be depleted onto dust. 

6. The constancy of IRX is consistent with the 
premise that neither grain destruction by sputter- 
ing nor cooling of the shocked gas have played a 
significant role in the evolution of the IR and X- 
ray emission. This scenario is consistent with the 
sputtering rates currently used in the literature as 
expressed in eq. (1). 

7. The presence of a secondary dust component, ra- 
diating at significantly higher temperatures than 
the silicate dust, suggests that the grain radii or 
IR emissivities of this component must be signif- 
icantly smaller than those of the silicates. Their 
sputtering lifetime could therefore be significantly 
shorter than that of the silicate grains and their 
evolution distinctly different from that of the sili- 
cates, especially around the 6100 d time interval. 
However, a significantly shorter grain lifetime con- 
tradicts the nearly constant or slightly decreasing 



flux ratio between the hotter secondary component 
and the silicate grain component. Hence, the na- 
ture of the secondary dust component remains a 
mystery. 

Five years of continued monitoring of SN1987A have 
lead us to conclude there is no current evidence for grain 
destruction. This conclusion stands in c o ntras t to the 
results of our earlier analysis iDwek et al.l (120081 ) , which 
were based only on two epochs of data at early times, and 
which we now see to be somewhat anomalous. Continu- 
ing observations will reveal evidence for grain destruction 
in the ER, and may elucidate the nature of the mysteri- 
ous secondary dust component. 
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